Coarsening Dynamics in a Two-dimensional XY model with Hamiltonian Dynamics 
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We investigate the coarsening dynamics in the two-dimensional Hamiltonian XY model on a 
square lattice, beginning with a random state with a specified potential energy and zero kinetic 
energy. Coarsening of the system proceeds via an increase in the kinetic energy and a decrease in 
the potential energy, with the total energy being conserved. We find that the coarsening dynamics 
exhibits a consistently superdiffusive growth of a characteristic length scale as Lit) ~ t 1 ^ with 
1/z > 1/2 (ranging from 0.54 to 0.57 for typical values of the energy in the coarsening region). Also, 
the number of point defects (vortices and antivortices) decreases as Nv{t) ~ t~^ v , with <f>v ranging 
between 1.0 and 1.1. On the other hand, the excess potential energy decays as AU ~ t~^ u , with a 
typical exponent of 4>u — 0.88, which shows deviations from the energy-scaling relation. The spin 
autocorrelation function exhibits a peculiar time dependence with non-power law behavior that can 
be fitted well by an exponential of logarithmic power in time. We argue that the conservation of 
the total Josephson (angular) momentum plays a crucial role for these novel features of coarsening 
in the Hamiltonian XY model. 

PACS numbers: 64.60.Ht,64.60.Cn,75.10.Hk,75.40.Gb 



I. INTRODUCTION 

In recent years, there have been research efforts toward 
investigating the statistical mechanical behavior of lattice 
spin systems by solving directly the Hamiltonian equa- 
tion of motion associated with the lattice spin system of 
interest. Especially, equilibrium phase transitions have 
been found to be associated with the existence of topo- 
logical change in the dynamics, as manifested in the dis- 
tributions of Lyapunov exponents^^iSiLSiS. 

In relation to these interests, we here investigate the 
coarsening dynamics of the XY model on a square lattice 
based on Hamiltonian dynamics. In conventional studies 
of phase ordering and coarsening dynamics in a statis- 
tical system, the system is quenched from the random 
disordered state to a low-temperature state below the 
transition temperatur e) 10 ! 11 . In computational studies, 
one employs Monte Carlo or Langevin dynamics meth- 
ods for simulating the coarsening processes. In typical 
situations, the average length scale L of ordered domains 
grows in time as a power law L(t) ~ t x l z , where the 
growth exponent 1/z depends on the dimension of the 
space and the dimension of the relevant order param- 
eter, in addition to the conserved or nonconserved na- 
ture of the latter in the relaxation dynamics^. Usually, 
characteristic topological defects, such as point vortices 
or domain walls, are generated in the initial disordered 
state, depending on the dimension of the order parame- 
ter, and the annihilation of these defects gives the main 
mechanism of coarsening and phase ordering in the sys- 
tem. The observed self-similarity of these coarsening sys- 
tems at different time instants is usually represented by 
the so-called dynamic scaling hypothesis of the equal-time 
spatial correlation function of the order paramete r^ 



In equilibrium, the XY model exhibits a Kosterlitz- 
Thouless (KT) transition at Tkt due to the unbinding of 
vortex-antivortex pairs 1 ^. Below Tkt, the system has a 
quasi-ordered phase that is characterized by a power-law 
decay of the order parameter correlation function for long 
distances. The critical exponent governing the power-law 
decay decreases continuously down to zero temperature: 
the system is critical at equilibrium for all temperatures 
below Tkt- 

After some efforts of various 

groupa 14 i 15 i 16 i 17 i 18 i 19 i 20 i 21 i 22 i 23 i 24 , as well as some con- 
troversies, it is now agreed that, in the phase ordering 
dynamics of the XY model via Monte Carlo simulations 
or Langevin dynamics methods, the growing length scale 
exhibits a logarithmic correction^ to diffusive growth as 
Lmc{P) ~ ^/(logi)] 1 / 5 . Here, the logarithmic correction 
can be attributed to a logarithmic divergence (in the 
system size) of the effective friction constant of a moving 
vortex in the dissipative dynamics with a non-conserving 
order parameter. In a related simulational work on 
the coarsening dynamics of superconducting juction 
arrays based on the dynamics of the resistively-shuntcd 
junction (RSJ) model, it was revealed that there is no 
logarithmic correction. This could be understood in 
terms of a finite friction constant of a moving vortex 
in the limit of large system size, which is due to the 
particular type of dissipative coupling in Josephson 
junction arrays^. 

One of the questions we would like to address in this 
work is whether we would find similar features as men- 
tioned above in the coarsening of the XY model based 
on Hamiltonian dynamics. Unlike the case of conven- 
tional coarsening based on dissipative dynamics, almost 
no works exist on coarsening via Hamiltonian dynamics, 
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except for the case of Ising-type systems. In the simpler 
case of the Ising-type </> model on a two-dimensional 
square lattice with nonconserved order parameter, coars- 
ening via Hamiltonian dynamics was investigated by 
Zheng 8 . He contends that the growth exponent is differ- 
ent from the usual diffusive growth of 1/z — 1/2 obtained 
in the traditional phase ordering. This result, however, 
was later disputed by others as a crossover effect-, which, 
we think, is still controversial because the asymptotic dy- 
namic features are not completely settled yet. Now, in 
the case of the XY model with Hamiltonian dynamics, 
it is of some interest to investigate whether any new fea- 
tures exist in the non-equilibrium coarsening processes. 
Due to the Hamiltonian dynamics nature of our system, 
the total energy is conserved, and we specify the initial 
states by their total energy divided by the system size, 
i.e., the per-site energy or the energy density. These ini- 
tial states with specified energies are prepared with spe- 
cial Monte Carlo algorithms to be explained in section 
II. In order to make an analogy to ordinary coarsening in 
dissipative dynamic relaxation more apparent, the initial 
rotational velocities of the rotors are taken to be zeroes. 

Since we begin with zero kinetic energy (zero rotational 
velocities for all spins) , the Hamiltonian dynamics of the 
system will generate kinetic energy taken from the po- 
tential energy. As time proceeds, we can expect that 
(for a per-site energy that is low enough to correspond 
to a low-temperature quasi-ordered phase), the system 
will evolve toward some equilibrium stationary dynamic 
state, which can be considered as the thermal equilib- 
rium state in the case of purely dissipative dynamics. In 
other words, the kinetic energy that is generated during 
the course of coarsening process acts as a thermal bath 
for the system. We investigate the time-dependent spin 
configuration for the system in terms of spatial ordering 
and relaxation of the vortex numbers, etc., in analogy to 
the conventional dissipative coarsening systems based on 
Monte Carlo or Langevin dynamics. 

We find that the equal-time spatial correlation func- 
tions approximately satisfy critical dynamic scaling 

C(r,t) = r-^g(r/L(t)), (1) 

with a spatial correlation exponent rj(T) (which is in- 
creasing in T) and a growing length scale L(t) that 
grows typically as L(t) ~ t 1 ^, where z is the dynamic 
exponent. In the late time region, we find that the 
length scale L(t) grows with an exponent 1/z that is 
larger than the diffusive exponent 1/2. Closely related 
to this growing length scale is the average separation be- 
tween vortices Ly(t), which can be derived from the de- 
cay of the number of vortices Nv(t), which behaves as 
Ny(t) — Nv(oo) ~ For unform distribution of vor- 

tices, we expect L(t) ~ L v (t) ~ (AV W)~ 1/2 ~ t~* v l 2 in 
the limit where Ay(oo) is negligible (which appears to be 
valid approximately in the low-energy region of E < 0.5). 
We see that the vortex density relaxations exhibit (ap- 
proximately) the expected behavior with the exponents 
4>v being larger than unity. 



We find also that the relaxation of the potential energy 
exhibits a power-law decay A/7 = U{t) — Uoo ~ t~^ u to- 
ward a long-time equilibrium value, with the exponents 
taking values approximately in the range 0.86 < 4>u < 
0.90 (in typical cases, (f>u ~ 0.87 - 0.88). This shows 
a considerable deviation from the energy-scaling rela- 
tion relating the excess energy and the growing length 
scaled. Another interesting feature of the coarsening dy- 
namics in the Hamiltonian XY model is that the spin 
autocorrelation exhibits a peculiar non-power-law be- 
havior, which can be reasonably well fitted by A(t) ~ 
Ao exp[— &(ln(i)) 7 ], with b ranging from 0.24 to 0.4 and 
the exponent 7 from 1.5 to 1.7. 

In order to understand the microscopic mechanism of 
these coarsening features, including the notable superdif- 
fusive growth exponent, that are quite distinct from the 
cases of purely dissipative dynamics, we note that one 
conserved quantity exists in this Hamiltonian model sys- 
tem. That is, the total angular momentum (or Josephson 
angular momentum) is conserved in the Hamiltonian dy- 
namics. Due to the conservation of angular momentum, 
propagating spin wave modes appeasS&Sii. We strongly 
believe that the existence of propagating spin wave modes 
due to the conservation of angular momentum is the main 
cause of the superdiffusive nature of the coarsening and 
of the other unique features in this system. In order 
to obtain an effective friction constant of a vortex, we 
could apply a collective variable approach. However, this 
method gave us a vanishing effective friction constant to 
the lowest order in the velocity of a vortex. Higher-order 
friction can be argued to increase the growth exponent 
1/z from the diffusive value of 1/2. Further work is nec- 
essary to reach a full understanding of the main features 
of the coarsening in the Hamiltonian XY model in two 
dimensions. 



II. THE HAMILTONAIN XY MODEL AND 
SIMULATION METHODS 



The potential energy function of the hard-spin XY 
model on a square lattice is given by 

[7 = 7^(1-008(^-^0), (2) 

<ij> 

where J is the interaction strength, 9i is the phase angle 
of the spin at site i and the sum is over nearest neighbor 
pairs. The corresponding Hamiltonian dynamic equation 
can be obtained by adding the kinetic term: 

ff=E^ + J E( 1 - cos (^-^)> (3) 

i <ij> 

where rrij = 9i represents the canonical momentum as- 
sociated with the angular variable 9i of planar spins at 
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site i. This model appears naturally in easy-plane ferro- 
magnets or superfluid helium with the conserved variable 
m corresponding to the z-component of the spins or the 
density of the superfluid respectively in the limit of neg- 
ligible themal noise (e.g., negligible phonon effect) 30 ' 31 . 

Previous works^ have shown that a Kosterlitz-Thouless 
(KT) transition occurs at the value of potential energy 
around Ukt — LO. Our computational scheme in this 
work performs Newtonian dynamic simulations based on 
the above Hamiltonian with initial states that take zero 
kinetic energy and specified potential energy values but 
are otherwise random. Specified potential energy values 
are chosen from the region that corresponds to below or 
near the KT transition energy (Ukt — LO). 

Unlike the case of dissipative systems with Langevin 
dynamics or Monte Carlo dynamics where dissipative re- 
laxation is incorporated in a natural way, here in our case 
of Hamiltonian dynamics, the total energy (sum of kinetic 
and potential energy) is conserved. However, since we be- 
gin the Hamiltonian dynamics with zero kinetic energy 
with a given specified potential energy, we expect that 
the system will evolve in such a way as to increase (on 
the average) the kinetic energy and decrease the potential 
energy. Eventually, in the long-time limit, the system will 
reach a certain kind of equilibrium dynamic state where 
the average potential energy and the kinetic energy do 
not vary anymore. Of course, in terms of phase-ordering 
dynamics in dissipative systems, a system quenched to a 
symmetry-broken phase will never reach complete equi- 
librium in a finite time period. 

Therefore, we can consider the kinetic part of the sys- 
tem acting as a kind of thermal heat bath that absorbs 
the " dissipated heat" coming from the decreasing poten- 
tial energy. We can, thus, investigate the time-dependent 
spin configuration for the system in terms of spatial or- 
dering and relaxation of vortex numbers, etc., in analogy 
to the Monte Carlo or Langevin relaxation dynamics. We 
generate the initial states with specified potential ener- 
gies in the following manner: To begin with, we take a 
random spin configuration. Then we employ a Monte 
Carlo annealing algorithm such that the spin configura- 
tions with potential energies that are closer to the target 
energy are preferentially accepted. In other words, we 
steer the initial state to some target spin configuration 
such that we obtain a state whose potential energy is 
equal to a specific energy; otherwise, the spin configura- 
tion is random. 

The Hamiltonian equation reads 



m. 



SU 



0i = m. 



-J^sin^-^), (4) 

3 

(5) 



where j denotes the nearest neighbors of site i and i takes 
values i = 1, • • ■ , TV 2 . For simplicity, we set the rotational 
inertia to be equal to unity and put J — 1. Here, we 
note that, in addition to the conservation of total energy, 
another conserved quantity, namely the total ( Josephson) 



angular momentum ^ 4 m t = 0i exists; This is due to 
the invariance of the equations under global translations 
in the variable 8i, that is, under 6i — > 0i + a. 

Equations 4 and 5 are numerically integrated in time 
using a second order velocity- Verlet algorithm 3 ^ with the 
time integration step of At = 0.01, which conserves the 
total energy to within a ratio of 10 up to the maximal 
time integration steps of around 10 6 . Periodic boundary 
conditions on both lattice directions are employed. The 
simulations are carried out on square lattices of dimen- 
sions up to N x TV = 1800 x 1800 with sample averages of 
over 20 to 40 different initial configurations. A parallel 
computation is employed on Linux clusters with from 16 
to 32 CPU's via a domain decomposition method with 
high parallel efficiency. 

The final results are obtained from averages over 10 
to 40 different random initial configurations. The main 
quantities of interest are as follows: 

(i) Equal time spatial correlation function of the order 
parameter, 



C(r,t) = i_ /^cos(0i(f) - 9 i+r (t)\ , (6) 

where < • • • > denotes an average over random initial 
configurations. 

(ii) Number of vortices Ny(t) at time t, including both 
positive and negative vortices. 

(iii) Time-dependent relaxation of the potential energy, 
U(t), 



= Ec 1 -™^-^))), ( 7 ) 



\<h3> 



(iv) Nonequilibrium spin autocorrelation function 



A(<)=^(x;«»(fc(o) -*(*)) 



(8) 



The results of simulations will be presented in section III. 



III. SIMULATION RESULTS AND 
DISCUSSIONS 



Figure 1 shows snapshots of spin configurations for the 
case of E = 0.4 at four different time stages of coarsening. 
We can see that coarsening occurs via decay of vortices 
and anti- vortices. Figure 2(a) shows the equal-time spa- 
tial correlation for E — 0.4. The simplest way to analyze 
the data is to attempt the so-called critical dynamic scal- 
ing 



C(r,t)=r-^g(r/L(t)), 



(9) 
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with the spatial correlation exponent rj(E) and a length 
scale L(t) growing typically as L(t) ~ t x l z . In our case, 
for a given time instant t, we determined L(t) in such 
a way that r-^C{r,t)\ {r=L{t)) = g(r / L{t))\ {r=L{t) ) = 
g(l) = 0.4. For a suitable choice of f]{E), we can obtain a 
good collapse of the rescaled correlation functions g(x), 
at least in the late-time region. Shown in Fig. 2(b) is 
such an attempt at critical dynamic scaling. We find 
that the equal-time spatial correlation functions exhibit 
a reasonable critical dynamic scaling, at least, in the late- 
time region. However, in the early-time regime, we can 
see that it is not possible to collapse the correlations with 
any reasonable values of 77. 



E 


(j>v {N = 1000) 


4>v (N = 1000) 


4>v(N = 1800) 


i> v (N = 1800) 


0.1 


1.001(6) 


0.870(15) 






0.2 


1.059(8) 


0.900(20) 


1.050(10) 


0.902(15) 


0.3 


1.085(10) 


0.879(15) 






0.4 


1.092(11) 


0.879(15) 


1.055(11) 


0.880(12) 


0.5 


1.060(11) 


0.867(15) 






0.6 


1.025 (15) 


0.867(16) 






0.7 






0.943(20) 


0.835(16) 


0.8 






0.950(20) 


0.890(18) 



TABLE II: Vortex number exponent <j>v and the potential 
energy exponent 4>u- 



E 


•q (N = 1000) 1/ 


z (N = 1000) 7i(7V = 1800) 1/ 


z (N = 1800) 


0.1 


0.046(10) 


0.520(12) 






0.2 


0.065(16) 


0.540(15) 


0.09(2) 


0.549(12) 


0.3 


0.104(16) 


0.566(10) 






0.4 


0.112(17) 


0.565(12) 


0.115(18) 


0.555(10) 


0.5 


0.120(18) 


0.549(13) 






0.6 


0.135 (19) 


0.549(13) 






0.7 






0.156(21) 


0.560(14) 


0.8 






0.170(22) 


0.565(14) 


1.0 






0.235(20) 


0.62(25) 



TABLE I: Correlation exponent r\ and the growth exponent 
1/z. 



In the late-time region, we find that the length scale 
L{t) grows with an exponent 1/z that is larger than the 
diffusive exponent 0.5. For example, in the case of E = 
0.4, we obtain 1/z ~ 0.555 ± 0.010, which is definitely 
larger than 1/2 (Fig. 2(c), Fig. 2(d), and Table 1). Note 
that the power law region shows up at t > tp « 50, which 
is a considerably late-time stage. 

Closely related to this growing length scale is the av- 
erage separation between vortices, Lv(t), which can be 
derived from the decay of the number of vortices, Nv(t). 
Shown in Fig. 3(a) are the number of vortices versus time, 
which behaves as Ny(t) — AV(oc) ~ t~^ v , where AV (00) 
denotes the limiting value of the vortex number at equi- 
librium. At low temperature, the equilibrium vortex den- 
sity becomes vanishingly small. We regard Nyi<x>) as a 
fitting parameter such that the relaxation of the excess 
vortex number obeys the best power-law behavior. We 
find that the relaxation of the vortex number exhibits an 
excellent power-law behavior except for the early-time re- 
gion where smooth plateau is seen before crossing over to 
the late-time power-law region. The exponent <f>v takes 
values in the range from 1.0 to 1.10, depending on the 
system size and the energy, except for the limit of higher 
per-site energy. For example, in the case of E = 0.4 with 
L = 1800, we get (j> v ~ 1.055. Table 2 and Figure 3(c) 
shows the independence of the (f>v 

Now, for a unform distribution of vortices, we can ex- 
pect the average separation between the vortices to sat- 



sfy (for negligible AV(oo)) 

Lit) ~ Lv(t) ~ (Avwr l/2 



t 



-4>v/2 



(10) 



t 1/z , we expect 1/z = <f> v /2. The 
d </V/2 obtained from simulations 



1/z is slightly larger than </>y/2. 
We find also that the relaxation of the potential energy 
hibits a power law decay AU = U(t) — U x ~ t~^ u 
toward a long-time equilibrium value, with the exponents 
taking values approximately in the range 0.86 < 4>u < 
0.90 as shown in Fig. 3(b) and Fig. 3(c). The value of 
this exponent is almost independent of the energy E. For 
the usual dissipative coarsening with a non-conservative 
order parameter in two dimensions, the excess energy 
relaxation can be related to the growing length scale Lit) 
as 



AU~L- 2 it) log(L(i)/r ), 



(11) 



where ro denotes the size of a vortex core2£. Interest- 
ingly, in our model system with Hamiltonian dynamics, 
we could not show confidently that the simulation re- 
sults are consistent with the above relation between the 
energy and the domain growth relation. Further work is 
necessary in this regard^. 

Another interesting feature of the coarsening dynamics 
in the Hamiltonian XY model is that the spin autocor- 
relation exhibits a peculiar non-power law behavior that 
can be reasonably well fitted by 



A(t) ~4)exp[-fc(ln(i)F], 



(12) 



with b ranging from 0.24 to 0.4 and the exponent 7 from 
1.5 to 1.7. One typical example is shown in Fig. 4 for 
the case of E — 0.4 with L = 1800, which shows an ex- 
cellent fit to the above functional form with b ~ 0.242 
and 7 ~ 1.67. Note that, in the limit of 7 = 1, Ait) re- 
duces to a power law of A{t) ~ Aot~ b . The fitted values 
of 7 (ranging from 1.5 to 1.7) implies that Ait) exhibits 
a considerable deviation from a power-law behavior. To 
the best of our knowledge, the Hamiltonian XY model 
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on a square lattice is the first system in which the au- 
tocorrelation function exhibits a functional form with an 
exponential of a logarithmic power. 

Since we have no theories available on the coarsening 
of the Hamiltonian XY model in two dimensions, we may 
resort to phenomenological approaches for the explana- 
tion of these results on the phase ordering of the Hamil- 
tonian XY" model, especially the relationship between the 
characteristic growth exponents and other exponents ob- 
tained from our simulations. 

To begin with, let us note that a strong contrast ex- 
ists between the above results of Hamiltonian coarsening 
and the conventional dissipative coarsening (via Monte 
Carlo or Langevin dynamics) in the XY model. The 
first thing to recall is that the growth law of the usual 
XY model with Monte Carlo dynamics (on nonconserved 
order parameter) follows Lit) ~ (t/ ln(i)) 1 / 2 . In con- 
trast, here, we observe in the Hamiltonian XY" model, a 
growth law with the exponent 1/z being greater than 
the exponent for the diffusive growth of 1/z = 1/2. 
Also, the vortex number density exhibits a relaxation 
N v - AV(oo) ~ t-* v with cj> v greater than 1.0. 

As we mentioned earlier, the system is expected to de- 
crease in its total potential energy via a corresponding 
increase of kinetic energy. This kinetic energy part is 
thought to consist of spin wave oscillations (superposed 
upon vortices or anti- vortices). We expect growth laws 
of length scales to be determined by the effective dynam- 
ics of vortices and anti-vortices and by the statistics of 
vortex-antivortex annihilation. 

Even though our model system is very simple in its 
equations, we find it not very easy to find an effective 
vortex dynamics for a given total energy of the system. 
One way may be to find a Langevin-type model system, 
which is expected to be equivalent to our microcanoni- 
cal Hamiltonian system in the thermodynamic limit, and 
then to investigate the coarsening dynamics of the equiv- 
alent model. 

Noting that the equation of motion can be derived from 
the Hamiltonian 



cos 



(13) 



with the Poisson bracket {9i, rrij} — g (where g — 1 here) 
and considering that the total sum of m^, i.e., J^. to^, is 
conserved, we can construct a Langevin equation for 9 
and m in the following way: We suppose, to begin with, 
that the fluctuating spin waves act as a thermal bath 
on the vortices, with the temperature determined by the 
equivalent equilibrium temperature corresponding to the 
total energy. Thus, we may suppose the effective coarsen- 
ing dynamics of the original Hamiltonian XY model sys- 
tem to be equivalent to the dissipative dynamics of the 
system under thermal noise associated with the equiva- 
lent temperature: 



5m. 
SH 

m * = -9T7T 



SH SH 
9— - X TT + W)> 



^V 2 ^- +&(£), 
om, 



where the Gaussian noises r\i and £i satisfy 
(Vi(t)Vj{t')) = 2Xk B T8 ij 8{t-t'), 

(&(*)&(0> = Wv\-i(t-0 



(14) 
(15) 

(16) 
(17) 



at temperature T, ks being the Boltzmann constant. 
Here, Sij represents the discrete delta function on a two- 
dimensional square lattice, and V 2 the corresponding dis- 
crete Lapacian reflecting the conservation of the total 
Josephson angular momentum, m. In the linear limit, 
ignoring the vortices, these set of equations reduce to the 
Nelson-Fisher model describing low-temperature spin dy- 
namics in one and two dimensions^. 

It is not a trivial matter to prove an equivalence be- 
tween the relaxation of the Hamiltonian XY model and 
that of the above Langevin equations with reversible 
mode-coupling. In order to be sure of the equivalence, 
we should assume that the spin dynamics is sufficiently 
chaotic so that the spin-wave fluctuations generate fully 
Gaussian random noise on the spin dynamics. Beginning 
with Eqs. 14 and 15, we may also obtain an effective dy- 
namics of a single- vortex excitation under external force. 

By assuming a slowly moving vortex, we attempted to 
apply a collective variable methodi^i&Si to Eqs. 14 and 
15 in order to obtain an effective equation of motion for a 
single vortex. We obtain from this method an interesting 
result that, due to the conservation of m, the effective 
equation for the slow-velocity limit of a vortex can be 
written as 



MV + CV = F + F n 



(18) 



with the effective mass M of a vortex behaving as M ~ 
ln(L), and the noise term F no i se being an uncorrelated 
Gaussian noise. £ represents the effective linear friction 
constant of the vortex. F is an external force, such as 
coulomb forces due to other vortices, etc. An interesting 
feature is that the dissipation term vanishes to the lowest 
order in the velocity of the vortex i.e., £ ~ + c V\ + ■ ■ ■ . 
Including the next order terms is a complicated task that 
we didn't actually carry out, but we can attempt a phe- 
nomenological argument on the coarsening by assuming 
a higher-order dissipation. Here, let us assume a dissi- 
pation term that is second order in the velocity of a vor- 
tex, and let us consider a vortex-antivortex pair with the 
coulomb force varying inversely in the mutual distance 
(in two dimensions). Then, a simple argument based on 
a force balance between the coulomb force and friction 
leads to (neglecting the inertial term in the limit of slow 
velocity and at low temperatures) 



cV 2 



K dt' r 



(19) 
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In this case, the time scale for annihilation of a vortex- 
antivortex pair with size R goes as t(R) ~ i? 3 / 2 . This, in 
turn, implies that the typical length scale between vor- 
tices will grow in time as L(t) ~ i 2 / 3 , which is definitely 
superdiffusive. This value of the growth exponent, 2/3, 
is a little larger than those exponents obtained from our 
simulations. This discrepancy might be attributable to 
the effect of vortex inertia, which diverges logarithmi- 
cally in the vortex size. Further work is necessary for a 
detailed quantitative understanding of the growth laws. 

We should note that, in our Hamiltonian dynamics 
simulations, the initial states are not completely random 
states due to the finite energy constraint, Therefore, the 
initial states already have some finite correlation length 
scale, i.e., separation L c between vortices. These length 
scales are of order of 3 to 5 lattice constants. Therefore, 
we expect that a finite time scale r c (corresponding to 
the time scale for a vortex to travel the distance L c ) ex- 
ists after which a scaling region emerges, which is well 
reproduced in our simulations. 



initial state that is specially tuned to have a given poten- 
tial energy (otherwise random) but with zero kinetic en- 
ergy develops into a late-time coarsening state, where the 
potential energy slowly decays as a power law with a com- 
pensating increase in the kinetic energy. We find that the 
coarsening dynamics exhibits a characteristic superdiffu- 
sive growth, with the exponents being consistently larger 
than the value for usual diffusive growth of 1/z = 1/2. 
Relaxation of various quantities cannot be understood in 
the usual framework of the dissipative dynamics of the 
Ginzburg-Landau Hamiltonian. We believe that these 
novel features of coarsening dynamics come from conser- 
vation of angular momentum, which, in turn, causes an 
effective nonlinear dissipation of individual vortices. 
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FIG. 1: Snapshots of the configuration of the XY spins in a 
square lattice of size L — 32, E = 0.4 and time t: (a) t = 0, 
(b) t = 0.01 • 2 7 , (c) t = 0.01 • 2 10 , and (d) t = 0.01 • 2 13 . 
Vortices and antivortices are denoted by solid squares and 
empty squares, respectively. 





FIG. 2: (a) Spatial behavior of the equal time XY spin cor- 
relation function at various time stages in square lattices of 
size N = 1800 at a per-site energy of E — 0.4. (b) Scaling 
collapse of the data in (a) with the appropriate scaling length 
Lit), (c) Scaling length versus time at various temperatures 
E = 0.4 and N = 1800, N = 1000. (d) Growth exponents and 
i](E) versus energy, manifesting superdiffusive growth. Error 
bars are of the order of the sizes of the symbols or smaller. 
Lines are only guides to the eyes. 




FIG. 3: Relaxation of (a) the vortex number N v at energy 
E = 0.40 and system size N = 1800 and 1000 and of (b) the 
excess potential energy AU at E = 0.4 in the case of size 
L = 1800. In (a), the power- law lines t~ 1Mi and i -1 ' 09 are 
also plotted whereas in (b) a power law line with ~ t~ ' 88 is 
shown, (c) shows the exponents <f>v and <j>u versus energy. 
Error bars are of the order of the sizes of the symbols. Lines 
are only guides to the eyes. 




FIG. 4: Spin autocorrelation function at energy E — 0.40 
and (b) the same in a double logarithmic plot. The fitted 
values are 6 ~ 0.242 and 7 ~ 1.67. 



